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The Monitor project: Rotation periods of low- mass stars in 
M50 
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ABSTRACT 

We report on the results of a time-series photometric survey of M50 (NGC 2323), a 
130 Myr open cluster, carried out using the CTIO 4m Blanco telescope and Mosaic- 
II detector as part of the Monitor project. Rotation periods were derived for 812 
candidate cluster members over the mass range 0.2 < M/Mq ^ l-l- The rotation period 
distributions show a clear mass-dependent morphology, statistically indistinguishable 
from those in NGC 2516 and M35 taken from the literature. Due to the availability 
of data from three observing runs separated by 10 and 1 month timescales, we are 
able to demonstrate clear evidence for evolution of the photometric amplitudes, and 
hence spot patterns, over the 10 month gap, although we are not able to constrain the 
timescales for these effects in detail due to limitations imposed by the large gaps in 
our sampling, preventing use of the phase information. 
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1 INTRODUCTION 

M50 is a populou s (~ 2100 stars brighter than V ~ 23; 



po p 

iKalirai et all 120031 ) open cluster of comparable age to the 



Pl eiades, at mo d erate distance. We adopt the parameters 
of iKalirai et"aLl l|2003h for the remainder of this work: a 
main sequence turn-off age of 130 Myr, distance 1000^*5 pc 
and redde ning E{B — V) = 0.22 mag (a n ave rage of val- 
ues f rom IClaria. Piatti fc Lapasseg Il998l and iHoag et al.l 
119671 : see IKalirai et all |2003| )." For comparison, the classi- 
cal m ain sequenc e turn-off a go for the Pleiades is 100 Myr 
(fMev net , Mermill iod fc Maed er 1993) . 

Despite having a very favourable angular size for CCD 
observations, there are only two previous CCD- based studies 
of the clust er population i n the literature, from IKalirai et al] 
(2003') and'Sh arma et "al] (|2006l ). obtaining similar values for 
the cluster parameters. We have undertaken a time-domain 
photometric survey of M50 with the joint aims of searching 
for eclipses due to stellar, brown dwarf and planetary com- 
panions orbiting low-mass cluster member stars, and to ob- 
tain a large sample of rotation periods to study the evolution 
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of stellar angular momentum in conjunction with previous 
measurements in the literature. 

In addition to studying angular momentum evo- 
lution, we can also begin to constrain the effects of 
other cluster parameters such as metallicity and en- 
vironment on the rotation period distributions, by 
comparing the M50 distribution with data for clus- 
ters of comparable age. Large samples of period mea- 
surements are available in NGC 2516 (~ 150 Myr; 
IJ effries. Thurston fc Hamblvl 120011 ) from llrwin et all 
(.2007bl ), cove ring 0.15 < M/Mp; < 0.7, and M35 (also 



150 Myr^ 



_ Barrado v Navascues. Delivannis fc Stauffeij 

I2OOII ) from iMeib om et al.l (|2008| ). covering masses down 
to ~ 0.6 Mq. Although there are relatively few rota- 
tion periods measured, especially for M-dwarfs, in the 
Pleiades, a large sample of v sin i measurem ents are avail- 
able c overing masses down to ^ 0.3 Mq JStauffer et al.l 



Il984l; IStauffer fc HartmannI Il987l ; ISoderblom et al 
Jon es. Fischer fc Staufferl 1 19961 : lOueloz et al.l 
Ter ndrup et all |2000|) . which can be compared to 



1993 



1998 



the 



rotation periods in a statistical sense. 

The rotational evolution of low-mass stars on the 
pre-main-sequence (PMS) is dominated by stellar contrac- 
tion, with angular momentum thought to be regulated 



© 0000 RAS 



2 J. Irwin et al. 



by processes relating to the star-disc int eraction (e.g. 
accretion-driven winds; Matt & Pudritz 2005, or "disc lock- 



ing^ |K6nig][l993; Collier Cameron, Camubcll fc Quaintrelll 
I1995I ), giving rise to a spread in rotation rates de- 
termined predomin an tly by the disc lifetime (e.g. 
Krishnamurthi et al.l Il997l: iBouvier. Forestini &: AUainl 



19971 : ISills. Pinsonneauh fc Terndrud 12000 "). As the stars 
arrive on the zero age main sequence (ZAMS), the stellar 
contraction ceases, and angular momentum losses via 
magnetised stellar winds dominate the subsequent evo- 
lution. Observations indicate that the rotation rates of 
solar-type stars betw een the age of the Hyades (~ 625 Myr; 
iPerrvman erall [l998l ) and the age of the Sun ar e well - 
described by u) PC t~^^^, the famous ISkumanichI (|l972l ) 
law (e.g. ISoderblomI 119831 ). This can be reproduced in a 
more theoretically-motivated framework from parametrised 
angular momentum loss laws (usually based on iKawalerl 
1 19881 ). 

Observations of the Pleiades and other young clusters 
at ~ 50 — 100 Myr represent a "snapshot" of the rotational 
evolution process at the point where solar-type stars have 
recently reached the ZAMS, and before significant angular 
momentum losses due to stellar winds have taken place. 
These indicate that although the evolution of the slowest 
rotators from this age to the age of the Hyades (and the 
Sun) can be reproduced by the lSkumanichi (1972) law, these 
clusters show a spread in rotation rates at a given mass, 
and a number of "ultrafast rotators" (e.g. ISoderblom et all 
1 19931 ). neither of which are seen in the Hyades and at older 
ages: here the rotation rate is typicall y found to follow a 
fairly well-defined function of mass (e.g. lRadick et al.lll987^ . 
In order to reproduce the ultrafast rotato rs on the ZAMS, 
most modellers modify the iKawalej 11988*) formalism to in- 
corporate saturation of the ang ular mo mentu m losses above 
a critical angular vel ocity Wsat (|Stauffer fc Hartmannlll987l : 
iBarnes &: Sofia|[l996l ). The saturation is further assumed to 
be mass-dependent, to account for the mass-dependent spin- 
down timescales observed on the early main sequence. 

There is mounting evidence in the literature that this 
picture of angular momentum evolution is still not sufficient 
to reproduce the observations on the ZAMS, in particular 
over the interval between ~ 50 — 500 Myr, if we assume 
the stars rotate as solid bodies. The observations indicate 
that such models produce a spin-down over this age range 
that is too rapid, particularly for th e slowest rotators in 
open clusters (e.g. Ilrwin et al.ll20"07bl). Several studies (e, 
iKrishnamurthi et all 1 19971 : lAllainlll998l : Ilrwin et all 12007 
have invoked core-envelope decoupling, where the radia- 
tive core and convective envelope of the star are allowed to 
have different rotation rates, as a means to produce a more 
shallow evolution, by coupling angular momentum from a 
rapidly-spinning core (which experiences little angular mo- 
mentum loss assuming the disc and wind couple predomi- 
nantly to the outer convective regions of the star) to the 
convective envelope on timescales of a few 100 Myr to pro- 
vide a "late time replenishment" of the surface rotation rate. 

Open clusters of 50 — 500 Myr age represent an ideal 
testing ground for these models. The rotational evolution is 
strongly mass-dependent, so in order to decouple the mass 
effect from the time-dependence, extremely large sample 
sizes are required, such that we can obtain good statistics 
over small bins in mass. M50 represents an ideal target for 



such a study, having a large population over a relatively 
small area of sky, which can be observed very efficiently us- 
ing the multiplex advantage afforded by a standard wide- 
field CCD mosaic. 

The remainder of the paper is structured as follows: the 
observations and data reduction are described in ^ and the 
colour magnitude diagram (CMD) of the cluster and candi- 
date membership selection are presented in ^ The method 
we use for obtaining ph otometric periods is summarised in 
El (see H rwm et al.ll2006l for a more detailed discussion). Our 
results are given in and ^summarises our conclusions. 



2 OBSERVATIONS AND DATA REDUCTION 

Photometric monitoring observations were o btained as part 
of the Monitor project l|Aigrain et al.l [20071 ). using the 4m 
Blanco telescope at Cerro Tololo Interamerican Observa- 
tory (CTIO), with the Mosaic-II imager. This instrument 
provides a field of view of ~ 36' x 36' (0.37 sq.deg), us- 
ing a mosaic of eight 2048 x 4096 pixel CCDs, at a scale of 
~ 0.27"/pix. A total of ~ 95 hours of photometric moni- 
toring was conducted, spread over three distinct observing 
runs: 6x full-nights in two three-night segments spanning 
2005 Feb 04-06 and 2005 Feb 14-16, 8 x 1/2-nights between 

2005 Dec 24 and 2006 Jan 06, and 4x full-nights between 

2006 Jan 28 and 2006 Feb 01. M50 was observed in a single 
telescope pointing centred on the cluster, for ~ 8 hours per 
night ('^ 4 hours per 1/2- night), in parallel with another 
field in the cluster NGC 2362, t he results from which were 
published in Ilrwin et al.l (|2008bh . Exposure times were 75 s 
in i-band, giving a cadence of ~ 6 minutes (composed of 2x 
75 s exposures plus 2x 100 s readout time, slewing between 
M50 and NGC 2362 during readout). We also obtained deep 
1^-band exposures (2 x 600 s, 450 s and 150 s) which were 
stacked and used to produce a colour-magnitude diagram of 
the cluster. 

For a full desc ription of ou r data reduction steps, the 
reader is referred to llrwin et al. I (l2007a). Brieflv. we used the 
pipeline for the INT wide-field survey (|lrwin fc Lewis|[200ll ) 
for 2-D instrumental signature removal (bias correction, fiat- 
fielding, defringing) and astrometric and photometric cali- 
bration. We then generated a master catalogue for each fil- 
ter by stacking 20 of the frames taken in the best conditions 
(seeing, sky brightness and transparency) and running the 
source detection software on the stacked image. The result- 
ing source positions were used to perform aperture photom- 
etry on all of the time-series images. We achieved a per data 
point photometric precision of ~ 2 — 4 mmag for the bright- 
est objects, with RMS scatter < 1 per cent for i ^ 19 (see 
Figure[T} . A signal-to-noise ratio of 5 (corresponding approx- 
imately to the detection limit for point sources on a single 
frame of the differential photometry) is reached at i 22.7. 

Our source detection software fiags any objects detected 
as having overlapping isophotes. This information is used, 
in conjunction with a morphological image classification flag 
also generated by the pipeline software (Irwin & Lewis 200l|) 
to allow us to identify non-stellar or blended objects in the 
time-series photometry. 

Photometric calibratio n of our data w as carried out us- 
ing regular observations of iLandoltl (|l992l ) equatorial stan- 
dard star fields in the usual way. 
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Figure 1. Plot of RMS scatter per data point (measured over the 
entire data-set) as a function of magnitude for the i-band observa- 
tions of a single field in M50, for all unblended objects with stel- 
lar morphological classifications. The diagonal dashed line shows 
the expected RMS from Poisson noise in the object, the diagonal 
dot-dashed line shows the RMS from sky noise in the photometric 
aperture, and the dotted line shows an additional 1.5 mmag con- 
tribution added in quadrature to account for systematic effects. 
The solid line shows the overall predicted RMS, combining these 
contributions. 

Light curves were extracted from the data for ~ 63 000 
objects, 42 000 of which had stellar morphological classifica- 
tions (~ 23 per cent of these are flagged as having overlap- 
ping isophotes by the source detection software, and thus 
may be blended), using our standard apertu re photometry 
techniques, described in llrwin et al.l (|2007al ). We fit a 2-D 
quadratic polynomial to the residuals in each frame (mea- 
sured for each object as the difference between its magnitude 
on the frame in question and the median calculated across 
all frames) as a function of position, for each of the 8 CCDs 
separately. Subsequent removal of this function accounts for 
effects such as varying differential atmospheric extinction 
across each frame. Over a single CCD, the spatially- varying 
part of the correction remains small, typically ~ 0.02 mag 
peak-to-peak. The reas ons for using this te chnique are dis- 
cussed in more detail in llrwin ilaLf lioOZi). 

For the production of deep CMDs, we stacked 20 i-band 
observations, taken in good seeing and photometric condi- 
tions, and all of the l/-band observations. The limiting mag- 
nitudes on these stacked images, measured as the approx- 
imate magnitude at which our catalogues are 50 per cent 
complete, were V ~ 24.4 and i ~ 23.6. 



Rotation periods of low-mass stars in M50 3 



AO FO KO MO M3 M5 




CD 



V- 

Figure 2. V versus V — I CMD of M50 from stacked images plot- 
ted as a "Hess diagram" (greyscale map of the number density 
of sources in 0.1 mag bins), for all objects with stellar morpho- 
logical classification. The cluster sequence is clearly visible on 
the right-hand side of the diagram. The boundaries of the region 
used to select photometric candidate members are shown by the 
dashed lines (all objects between the dashed lines were selected). 
The reddening vector for Ay = 1.0 is shown at the right-hand 
si de of the diagram. T he mass scale is from the NextGen models 
of iBaraffe et ah I lll998l) . interpolated to 130 Myr, using our em- 
pirical isochrone to convert the V magnitudes to / magnitudes, 
and subsequently obtaining the masses from these, due to known 
pro blems with the V m agnitudes from the models (see, for exam- 
ple, H^afeeTal] Hasl) . The error bars at the left-hand side of 
the plot indicate the typical photometric error for an object on 
the cluster sequence. 



3 SELECTION OF CANDIDATE LOW-MASS 
MEMBERS 

3.1 The V versus V -I CMD 

Our CMD of M50 is shown in Figure [H The V and i mea- 
surements were converted to the standard Johnson- Cousins 
photometric system using colour equations derived from our 
standard star observations: 

(V-/) = (Kcd - iced) / 0.899 (1) 
V = Kcd + 0.005 {V - /) (2) 



/ = iced - 0.096 {V ~ I) (3) 

Candidate cluster members were selected by defining 
an empirical cluster sequence 'by eye' to follow the clearly- 
visible cluster single-star sequence. The cuts were defined by 
moving this line along a vector perpendicular to the cluster 
sequence, by amounts k — a{V — I) and k + a{V — I) as 
measured along this vector, where a{V — I) is the photo- 
metric error in the V — I colour. The values of k used were 
—0.125 mag for the lower line and 0.25 mag for the up- 
per line on the diagram, making the brighter region wider 
to avoid rejecting binary and multiple systems, which are 
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Figure 3. Contamination estimated from Galactic models, mea- 
sured as the ratio of the calculated number of objects in each 
magnitude bin from the models, to the number of objects de- 
tected and classified as candidate cluster members in that magni- 
tude bin. Note that bins with contamination estimates > 100 per 
cent (where there were more objects in that bin from the Galactic 
model than were actually observed) have been truncated to 100 
per cent. 



overluminous for their colour compared to single stars. 4249 
candidate photometric members were selected, over the full 
V magnitude range from V — 15.5 to 26, but the cluster 
sequence becomes difhcult to distinguish from the field pop- 
ulation for V'>24 {M<0.2 Mq). 

In this and subsequent sections, we make use of mass 
estimates for the cluster members. These were derived 
from the 7-band ab solute magnitudes and the models of 
iBaraffe et al] (|l998l '). interpolated to 130 Myr age. We did 
not use the V — I colour due to known problems with the V - 
band magnitudes from the models (see lBaraffe et aLlll998h . 



3.2 Contamination 

Since the cluster sequence is not well-separated from the 
field in the CMD, it is important to estimate the level of 
field star contamination in the sample of candidate cluster 
members. In this work, we have used Galactic models to ob- 
tain a simple first estimate of the contamination level. This 
will be refined using spectroscopic follow-up observations in 

a later publication. 

The Besangon Galactic models (|Robin et al.1120031 ) were 
used to to generate a simulated catalogue of objects pass- 
ing our selection criteria at the Galactic coordinates of 
M50 {I = 221.7°, b = -1.3°), covering the total FoV of 
~ 0.35 sq.deg (including gaps between detectors). The same 
selection process as above for the cluster members was ap- 
plied to this catalogue to find the contaminant objects. A 
total of 1970 simulated objects passed these membership se- 
lection criteria, giving an overall contamination level of ~ 47 
per cent after correcting for bins where the number of ob- 
jects predicted by the models exceeded the number actually 
observed (we simply assumed 100 per cent field contamina- 
tion in these bins). Figure |3] shows the contamination as a 
function of V magnitude. Note that this figure should be 
treated with extreme caution due to the need to use Galac- 
tic models, and especially given the overestimation of the 
numbers of observed objects by the models. 



4 PERIOD DETECTION 
4.1 Method 

The method we use for detection o f periodic variables was 
described in detail in llrwin et al.l (|2006l ). It uses least- 
squares fitting of sine curves to the time series m{t) (in 
magnitudes) for all candidate cluster members, using the 
form: 



m{t) — mdc + a. siniujt -f- 1 



(4) 



where mac (the DC light curve level), a (the amplitude) and 
<j) (the phase) are free parameters at each value of lu over an 
equally-spaced grid of frequencies, corresponding to periods 
from 0.05 — 50 days for the present data-set. 

For M50, we modified this to fit separate coefficients 
rridc, 01 and (j> for each of the three observing runs (a total 
of 9 coefficients at each value of u)) to allow for changes in 
the spot patterns that give rise to the observed photometric 
modulations. The period was required to remain the same 
over all the observing runs, which is expected since it should 
represent the underlying rotation period of the star, and 
this will change by a negligible amount due to rotational 
evolution over only 1 yr. This procedure is necessary for 
such a long data-set since significant evolution of the spot 
coverage of our targets is expected over timescales shorter 
than 1 yr. 

Periodic variable light curves were selected as before by 
evaluating the change in reduced x^: 



Ax. 



'X.f Xi^, smooth ^ 0.4 



(5) 



where Xv is the reduced x of the original light curve with 
respect to a constant model, and X^, smooth is the reduced 
X^ of the light curve with the smoothed, phase-folded ver- 
sion subtracted. Again, xS.smooth calculated for the three 
observing runs separately, and then summed, to account for 
any evolution in amplitude and/or phase. Since fitting for 
separate rridc coefficients in each observing run renders the 
period search insensitive to very long-period modulations 
appearing as simple "DC offsets" between the runs, we eval- 
uated xS for a three constant model, where we allowed for 
a different constant light curve level in each observing run. 

The threshold of Ax^ > 0.4 was used for the M34 data 
and appears to work well here too, checked by examining all 
of the light curves for two of the detectors, chosen randomly. 
A total of 1700 objects were selected by this automated part 
of the procedure. These light curves were then examined by 
eye, to define the final sample of periodic variables. A total of 
812 light curves passed this final stage, where the remainder 
had spurious variability (caused by systematic effects) or 
were too ambiguous to be included. 



4.2 Simulations 

Si mulations w e re per formed following the method detailed 
in llrwin et al.l (|2006l ). injecting simulated signals of 2 per 
cent amplitude and periods chosen following a uniform dis- 
tribution on logj^Q period from 0.1 to 20 days, into light 
curves covering a uniform distribution in mass, from 1.0 to 
0.1 M0. A total of 2111 objects were simulated. The phase 
of the modulations was randomised for each observing run in 
order to provide a more realistic evaluation of the sensitivity 
of our method across changes in spot patterns. 
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Figure 5. Detected period as a function of actual (input) pe- 
riod for our simulations. Objects plotted with crosses had frac- 
tional period error < 10 per cent, open circles > 10 per cent. 
The straight lines represent equal input and output periods, and 
factors of 2, 3, 1/2 and 1/3. The curved lines are the loci of the 
±1 day^^ aliases resulting from gaps during the day. The major- 
ity of the points fall on (or close to) the line of equal periods. 

The results of the simulations are shown in Figure |4] as 
diagrams of completeness, reliability and contamination as 
a function of period and stellar mass. Broadly, our period 
detections are close to 100 per cent complete for these ampli- 
tudes from 1.0 Mq down to 0.3 Mq, with remarkably little 
period dependence. At the very lowest masses (particularly, 
in the 0.2 < M/Mq < 0.3 bin of the diagram), the complete- 
ness drops substantially, and the reliability of the detected 
periods deteriorates slightly, due to the increased noise level 
in the light curves. For AI < 0.2 Mq, the completeness drops 
essentially to zero for 2% modulations. 

Figure [S] shows a comparison of the detected periods 
with real periods for our simulated objects, indicating good 
reliability of period recovery, with only a small amount of 
aliasing, especially when compar ed to some of our previous 
data-sets (e.g. M34; llrwin et al.i i2006). presumably due to 
the extended time-coverage of the M5Q observations. 



5 RESULTS 

The locations of our detected periodic variable candidate 
cluster members on a V versus V — I CMD of M5Q are 
shown in Figure [6] The diagram indicates that the majority 
of the detections lie on the single-star cluster sequence, as 
would be expected for rotation in cluster stars as opposed 
to, say, eclipsing binaries. 

The properties of all our rotation candidates are listed 
in Table [T] 

5.1 M50 rotation periods 

Plots of period as a function of V — I colour and mass for 
the photometrically-selected candidate cluster members are 
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Figure 6. Magnified V versus V—I CMD of M50, for objects with 
stellar morphological classification, as Figure |2] showing all 812 
candidate cluster members with detected periods (black points). 
The dashed lines show the cuts used to select candidate cluster 
members (see ^3.1|l . 

shown in Figure[7l Below ~ 0.7 Mq (or MO), these diagrams 
reveal a correlation between stellar mass (or spectral type) 
and the longest rotation period seen at that mass, with a 
clear lack of slow rotators at very low masses. This trend 
is also followed by the majority of the rotators in this mass 
range, with only a tail of faster rotators to ~ 0.2 day periods, 
and very few objects rotating faster than this. This is very 
similar to wh at we found in the earlier NGC 2516 a nd NGC 
2547 studies l|lrwin et al.ll2007bl : llr^in et al.ll2008al ). 

Above ~ 0.7 Mq, the M50 data indicate an inverse 
trend, of decreasing rotation period (faster rotation) as a 
function of increasing mass. The slope of this relation is 
much shallower than the one bel ow 0.7 Mq. The existence 
of such a relation was noted by iHartman et al.l (|2008l ) in 
M37 550 Myr) , and is also clearly evident in the v sin i 
data for NGC 2516 and the Pleia des, and re cently published 
rotation period data in M35 (Meibom et al.i r2Q08). as shown 
in Figure [TTI (left panel), which will be discussed in ^5.3\ 

These morphological features do not appear to be a re- 
sult of sample biases. In particular, the simulations show 
that the survey is sensitive to much shorter periods than 
0.2 day, and the upper limit in detectable periods is not 
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Figure 4. Results of the simulations for 0.02 mag amplitude expressed as percentages, plotted as a function of mass (left) and period 
(right). The simulated region covered 0.1 < M/Mq < 1.0 in order to be consistent with the M50 sample. Top panels: completeness as a 
function of real (input) period. Centre panels: Reliability of period determination, plotted as the fraction of objects with a given true 
period, detected with the correct period (defined as differing by < 20 per cent from the true period). Bottom panels: Contamination, 
plotted as the fraction of objects with a given detected period, having a true period differing by > 20 per cent from the detected value. 



Identifier RA Dec V o"(V) / <'"(^) P cti cr{cti) ci'z 0(012) «3 ""(as) 

J2000 J2000 (mag) (mag) (mag) (mag) (days) (mag) (mag) (mag) (mag) (mag) (mag) (M(r 



M50-1-98 


07 01 35.53 


-08 33 10.6 


17.082 


0.001 


15.911 


0.001 


12.807 


0.0091 


0.0003 


0.0035 


0.0002 


0.0115 


0.0011 


0.8 


M50-1-413 


07 01 38.56 


-08 34 35.7 


21.993 


0.015 


18.950 


0.007 


0.741 


0.0078 


0.0010 


0.0209 


0.0011 


0.0229 


0.0012 


0.3 


M50-1-651 


07 01 40.52 


-08 33 28.3 


16.521 


0.001 


15.423 


0.001 


6.811 


0.0055 


0.0002 


0.0063 


0.0002 


0.0024 


0.0004 


0.8 


M50-1-979 


07 01 43.77 


-08 32 58.5 


21.901 


0.014 


18.710 


0.006 


0.865 


0.0189 


0.0009 


0.0108 


0.0009 


0.0104 


0.0014 


0.4 


M50-1-1130 


07 01 45.15 


-08 34 11.0 


21.330 


0.009 


18.538 


0.006 


0.924 


0.0138 


0.0014 


0.0173 


0.0008 


0.0292 


0.0013 


0.4 



Table 1. Properties of our 812 rotation candidates, including V and /-band magnitudes and uncertainties from our CCD photometry, 
the period P, ^-ban d amplitudes ai, 02 , «3 and uncertainties (magnitudes, in the instrumental bandpass), interpolated mass and radius 
(from the models of lBaraffe et al.lll998l . derived using the / magnitudes). In the table, aix) denotes the uncertainty in quantity x. Note 
that these uncertainties do not incorporate systematic errors (e.g. in the zero point calibration), that dominate at the bright end. Our 
identifiers are formed using a simple scheme of the cluster name, CCD number and a running count of stars in each CCD, concatenated 
with dashes. The full table is available in the electronic edition. Machine readable copies of the data tables from all the Monitor rotation 
period publications are also available at |http : //www ■ ast ■ cam ■ ac . uk/research/monlt or /rotation/ [ 



mass-dependent, so this cannot explain the morphology of 
the upper envelope of rotation periods in Figure [7] More- 
over, the presence of these features in multiple rotation pe- 
riod data-sets with differing selection biases, and in v sin i 
data, strongly indicates that they are real features of the 
underlying rotation rate distribution of low-mass stars. 

As a further check, we can examine the distribution of 
rotation periods as a function of amplitude. This is shown in 
Figure |8] split into two mass bins. There does not appear to 
be any clear sensitivity bias, with very low amplitudes being 



detected over the full range of period, for the M < 0.5 Mq 
bin. 

The higher-mass bin (0.5 < A'I/Mq < 1.1) appears to 
show a correlation between amplitude and rotation period, 
in the sense that longer-period objects show smaller am- 
plitudes. This is confirmed by applying a non-parametric 
Spearman rank correlation analysis, which gave Vs = —0.25, 
corresponding to a probability of 3 x 10~® that the quan- 
tities are uncorrelated given the sample size. This trend is 
also present in the lower-mass bin, but at a lower signifi- 
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Figure 7. Plots of rotation period as a function of dereddened 
V — I colour (top), and mass (bottom) for M50, deriving the 
masses using the 13 Myr NextGen mass-magnitude relations of 
iBaraffe et"^ l ll998l) and our measured 7-band magnitudes. In the 
lower diagram, the greyscales show the completeness for 0.02 mag 
periodic variations from the simulations. 



cance level (mostly due to the smaller sample size), with 
Ts = —0.22, and a probability of 4 x 10"''. 

This correlation could be the result of a bias against 
detection of very small-amplitude modulations at short pe- 
riods, but this is in the opposite sense to the expected trend: 
it is easier to detect small modulations at short periods since 
the slope of the light curve within each night is larger, which 
contributes significantly to our confidence that the variation 
is real. We speculate that there are two likely causes of the 
trend we have observed: either this is a result of the well- 
known r elation betwe en rotation rates and stellar activity 
(e.g. see iGiidell [20o3 for a review), in the sense that fast 
rotators are more active, and therefore are more spotted, 
giving rise to a larger photometric amplitude; or, the con- 
tamination of our sample by field stars, which are in general 
less active and more slowly rotating than cluster stars, leads 
to a population of contaminant objects at the long-period 
end giving rise to the apparent correlation. The upper panel 
of Figure[8]already contains a hint that the latter may be the 
case: the distribution of amplitudes for periods > 10 days 
appears to be systematically skewed toward smaller ampli- 
tudes than that for periods < 10 days. Moreover, we expect 
there to be very few cluster stars rotating more slowly than 
10 days a t this age from rota tional evo lution mo dels (e.g. see 
llrwin et al. 2007bi : llrwin et al. 2008ai : llrwin et al.„2008bi ). 
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Figure 8. Plot of amplitude as a function of period for M50 
in two mass bins: 0.5 < M/M© < 1.1 (top) and M < 0.5 Mq 
(bottom). 



5.2 Radial distribution of periodic variables 

In order to further examine the effect of field contamination, 
we show in Figure [9] the radial distribution of the detected 
periodic variable objects, and in Figure [TO] the rotation pe- 
riod distribution as a function of mass in three radial bins, 
measuring the radius r with respect to the approximate po- 
sition of the cluster centre. We adopt a positio n of R.A. 
7''02™47.4=, Dec. -8°20'43" (|Sharma et al.ll2006l : see their 
table 5). These authors derive a cluster core radius of 6.5' 
and a total extent of 17', which are broadly consistent with 
Figure |9l so we expect field stars to dominate our largest- 
radius bin in Figure 1101 While this is clearly not the case, 
given the presence of the same well-defined morphology in 
rotation period as a function of mass in all three bins (pre- 
sumably due to the cluster population), the fraction of ob- 
jects with periods > 10 days does appear to increase with 
r, becoming especially clear in the bottom panel of the fig- 
ure. This suggests that these objects may indeed be slowly- 
rotating field stars contaminating the sample. 

The existence of a significant population of cluster stars 
in the lower panel of Figure [TO] is not necessarily inconsis- 
tent with the conclusions of Sharma et al.. (,2006, 1 . Firstly, by 
requiring the objects to have measured rotation periods, we 
strongly bias the sample in favour of cluster stars, since these 
are more active. Secondly, many clusters are found to have 
a larger appare nt exten t at low masses (e. g. lLarsonlll982l : 
iMcNamara fcSekiguchii 1 19861: ISagar et all Il988l ). in some 
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cases even havi ng a "halo" population of low-m ass stars (e.g. 
h and x Persei: ISchildlll967l : ICurrie et al.ll2007h . Our rotation 
period sa mple probes ~ 3 mag deeper in l/-band than the 
survey bv lSharma et al.) (|2006l ). and is therefore dominated 
by lower-mass stars. The radial distribution shown in the 
lowest-mass bin of their Figure 4 supports this hypothesis, 
showing a larger cluster population at large radii than their 
higher-mass bins. Due to the limited angular coverage of the 
present survey, it is difficult to confirm this from Figure (5] 



Figure 9. Radial distribution for objects with detected periods. 
The counts in each annulus have been corrected for the area of sky 
covered by the CCD mosaic by summing the area of all the pixels 
in falling within that annulus. This correction is only approximate 
for the outermost regions (e.g. r>0.35°) since it does not take 
into account incompleteness in the detections close to the edges 
of the detector, so the counts in the 0.35° < r < 0.4° bin in 
particular may be unreliable. The apparent deficit of objects in 
the innermost (r < 0.05°) bin is most likely due to a small error 
in the position of the cluster centre. 
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Figure 10. Rotation period as a function of mass, as Figure [S] 
plotted in throe radial bins of approximately equal area on the 
sky: r < 11' (top), 11' < r < 16' (centre), and r > 16' (bottom). 
The latter bin extends to the edge of the ~ 36' X 36' field of view 
of our observations. 



5.3 Comparison with NGC 2516 and M35 

By comparing the rotation period distributions in different 
clusters of the same age, we can begin to constrain the ef- 
fect of cluster environment, and other parameters such as 
metallicity, on the distribution of rotation periods. In the 
case of M50, suitable rotation period s amples are available 
in the literature for N GC 2516 jlrwin et al ., 200 7bi) and M35 
jMeibom et ahllioOSl ). with a reasonable overlap in the mass 
range covered. The metallicity of M50 does not yet ap- 
pear to have been reported in the literature, but estimates 
are available for NGC 2516 indicating that i t is near-solar 
([Fe/H] = -0.05 ± 0.14: lTerndrup eraLll2002l ). and for M35 
indicating a sub-solar metallicity ([Fe/H] = —0.21 ± 0.10; 
iBarrado v Navascues et al.l I2OOII ) . Figure [TT] shows a com- 
parison of the rotation periods measured in these two clus- 
ters with the M50 distribution from the present study. 

Examining the left-hand panels of Figure [TT] first, the 
morphology displayed in all the clusters is clearly very sim- 
ilar over the mass ranges in common. The Pleiades t;sini 
data appear to show a slightly different "turn-over" point 
for the slowest rotators, at ~ 0.9 Mq rather than 0.7 M0 
for the other clusters, but this may be the result of the sini 
ambiguity and the difficulty of detecting the slowest rotating 
objects via this method. Therefore, qualitatively, the distri- 
butions of rotation period versus mass for all the clusters 
appear to be very similar. 

In order to allow a quantitative comparison to be made, 
we have compared the distributions of rotation periods in 
three mass bins, chosen to resolve morphological changes in 
the diagram while concentrating on the regions where the 
various samples overlap sufficiently to allow a meaningful 
comparison. These are shown in the right-hand panels of 
Figure [TT] In all three cases, the distributions are extremely 
similar, generally differing by less than their combined re- 
spective Poisson counting uncertainties. This was confirmed 
by applying non-parametric two-sided Kolmogorov-Smirnov 
tests to the distributions, the results of which are given in 
the figure. In all cases the distributions are statistically in- 
distinguishable. 

The present samples in these three clusters therefore do 
not show evidence for dependence of rotation on cluster en- 
vironment, or a "third parameter" such as metallicity. This 
is of course not conclusive evidence that such a dependence 
does not exist, especially since the environment and other 
parameters are relatively similar for the three clusters con- 
sidered in this work. In order to conduct a more stringent 
test, additional rotation period observations in a substan- 
tially different environment at similar age will be needed. 
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rotation p eriod distribution wi th rotation period measurements in the Pleiades 
jVan Leeuwen. Alphenaar fc M cvs 1987: IStauifer et al.lll 987: Magnitskii IQsjl: iProsser et al.lll993al : iProsser et al.|[l993bl: iProsser et al I 



19951: iKrish namurt hi et ahl 1998i: iTerndrup et al.ll99g|:IScholz fc Ei slofFcl 200 j), v sin i measurements in the Pleiades llStauffer et al.lll984l : 
Stauffer fc Hartmann| ]l987l : ISoderblom et al.lll993l: Ijones. Fischer'fc Stauffcrl fT99^ : lOueIoz et ai]|l998l : iTerndrup et aI.ll2000^ ■ and rota- 



tion periods in NGC 2516 and M35. Left: rotation period as a function of ma ss for all five samples. Masses were computed using the 
measured 7-band magnitudes, interpolating the models of Bar affe et al. to the appropriate metallicity and age for each cluster, 

using the estimates of these parameters from the literature quoted in sections [Tl and [5.31 (for M50, we simply assumed solar metallicity 
given the lack of an estimate). Right: comparison of histograms of rotation period binned in mass to those regions where the cluster 
samples have sufficient overlap for a meaningful comparison. The grey shaded regions in the left panels indicate the mass ranges used 
for the histograms on the right. Note that these distributions have not been corrected for contamination by field stars due to the lack 
of follow-up data. Nevertheless, the rotation period distributions in all these clusters are qualitatively very similar, and appear to be 
statistically indistinguishable given the present sample sizes. 



5.4 Changes in amplitude: evolution of the spot 
patterns 

The availability of three epochs of data corresponding to 
the three different observing runs allows us to constrain the 
evolution of the spot patterns giving rise to the photometric 
modulations, by examining the evolution of the amplitude 
between the runs. Moreover, by comparing results in mul- 
tiple clusters (e.g. our NGC 2362 data-set, which has the 
same sampling) , we can examine the dependence of the spot 
evolution rate on cluster age. 

Figure [12] shows histograms of the ratio of the ampli- 
tudes between pairs of observing runs. As expected, over 
the interval between the first two runs (~ 10 months), com- 
pared to that between the second two runs (~ 1 month), 
the standard deviation of the amplitude ratio is larger, by 



a factor of ~ 1.8, indicating that a larger fraction of the 
amplitudes evolved over the longer time-span between the 
first two observing runs. This clearly demonstrates that the 
spot patterns on our target stars must have evolved over 
the 10 month gap. The evolution in amplitude is, nonethe- 
less, relatively modest, which is not surprising since it is 
reasonable to expect stars of a given activity level to main- 
tain comparable levels of spot coverage over long time-scales, 
with the phase of the modulations evolving according to the 
short time-scale evolution of individual spots or spot groups. 

Unfortunately, for the sampling strategy we have used, 
we cannot use the phase information to examine this evolu- 
tion. This arises because the number of cycles of modulation 
over the ~ 10 month gap between the first and second ob- 
serving runs is large, and unknown. Therefore, small errors 
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Figure 12. Histograms of the relative amplitudes from the sine 
curve fitting procedure described in i|4. 1 1 for the second and first 
observing runs (upper panel), and the third and second observing 
runs (lower panel). The mid-points of these were separated by 
~ 10 months and 1 month respectively. The open histogram shows 
the entire sample, and the solid histogram only those objects with 
amplitudes > 1%. The standard deviations of the distributions, 
in logjQ units, for the open (solid) histograms are 0.25 (0.17) for 
the upper panel, and 0.14 (0.10) for the lower panel. 




02 - 03 (rad) 

Figure 13. Histograms of the relative phase from the sine curve 
fitting procedure described in i|4.1l for the second and first observ- 
ing runs (upper panel), and the third and second observing runs 
(lower panel). The upper distribution is consistent with being uni- 
form, indicating that the phase is essentially randomised over the 
10 month time-span between the first and second observing runs. 
The standard deviation of the lower distribution is ~ 1.0 rad. 



in the rotation period, accumulated over the many cycles oc- 
curring in the 10 month gap, can lead to large phase shifts 
between the two observing runs. For example, even a one 
per cent error in a rotation period of 10 days gives rise to a 
cumulative uncertainty of ~ 0.3 cycle after 10 months, and 
this effect is worse for shorter periods. Figure [T51 illustrates 
this effect, showing that over the gap between the first and 
second observing runs, the phase is essentially randomised 
for 100% of the sample. 

Therefore, with the present data-set, it is difficult to fur- 
ther constrain the spot lifetimes giving rise to the observed 
evolution. 



6 CONCLUSIONS 

We have reported on results of an i-band photometric sur- 
vey of M50, covering ~ 0.4 sq.deg of the cluster. Selec- 
tion of candidate members in a. V versus V — I colour- 
magnitude diagram using an empirical fit to the cluster se- 
quence found 4249 candidate members, over a V magni- 
tude range of 15.5 < V < 26 (covering masses in the range 



0.1 < M/Mq < 1.1). The likely field contamination level was 
estimated using a simulated cat alogue of field obj ects from 
the Besangon Galactic models l|Robin et al.l 12003 ). finding 
an overall contamination level of ~ 47 per cent, implying 
that there are ~ 2300 real cluster members over this mass 
range in our field-of-view. 

We derived light curves for ~ 63 000 objects in the M50 
field, achieving a precision of < 1 per cent per data point 
over 15 < i < 19. The light curves of our candidate cluster 
members were searched for periodic modulations, presum- 
ably due to stellar rotation, giving 812 detections over the 
mass range 0.2 < M/Mq < 1.1. 

The rotation period distribution as a function of mass 
was found to show a clear mass-dependent morphology, sta- 
tistically indistinguishable from the distributions in the liter- 
ature for M35 and NGC 2516 once the different mass ranges 
probed by the surveys are taken into account. Thus the M50 
results do not yet provide any indication of a dependence of 
rotation rates on cluster environment, or a "third parame- 
ter" such as metallicity. 

Finally, we demonstrated evidence that the photometric 
amplitudes of a significant fraction of our targets do indeed 
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evolve over timescales of 10 months, as expected (for ex- 
ample, by analogy with the Sun). 
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